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Chapter 8 Volumes of parallelograms 


In the present short chapter we are going to discuss the elementary geometrical objects 
which we call parallelograms. These are going to be contained in some ambient Euclidean 
space, but each will have its own innate dimension. We shall usually denote the ambient space 
as R%, and we use its usual inner product e and norm. In discussing these parallelograms, 
we shall always for simplicity assume one vertex is the origin 0 € R%. This restriction is of 
course easily overcome in applications by a simple translation of RY. 






































DEFINITION. Let 7,...,2, be arbitrary points in R* (notice that these subscripts do 
not here indicate coordinates). We then form the set of all linear combinations of these points 


of a certain kind: 


i=1 





We call this set an n-dimensional parallelogram (with one “vertex” 0). We also refer to the 
vectors £1,...,% as the edges of P. O 


e For n = 1 this “parallelogram” is of course just the line segment |0, x1]. a ae 


X4 O 


e For n = 2 we obtain a “true” parallelogram. 


e For n = 3 the word usually employed is parallelepiped. ‘ ie 
However, it seems to be a good idea to pick one word *j 
to be used for all n, so we actually call this a 
3-dimensional “parallelogram.” 
wg 
We are not assuming that the edge vectors 71,..., 2, are linearly independent, so it could 


happen that a 3-dimensional parallelogram could lie in a plane and itself be a 2-dimensional 
parallelogram, for instance. 











PROBLEM 8-1. Consider the “parallelogram” in R* with “edges” equal to the three 
points 7, 7, 7? — 27. Draw a sketch of it and conclude that it is actually a six-sided figure 
in the x — y plane. 





A. Volumes in dimensions 1, 2, 3 
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We are going to work out a definition of the n-dimensional volume of an n-dimensional 
parallelogram P, and we shall denote this as vol(P); if we wish to designate the dimension, we 
write vol,(P). In this section we are going to recall some formulas we have already obtained, 
and cast them in a somewhat different format. 

e Dimension 1. We have the line segment [0,7] C R%, and we simply say its 1-dimensional 
volume is its length, i.e. the norm of x1: 





vol([0,21]) = ||x1|. 


e Dimension 3. We skip ahead to this dimension because of our already having obtained 
the pertinent formula in Section 7C. Namely, suppose P C R°, so we are dealing at first with 
the dimension of P equal to the dimension of the ambient space. Then if we write 71, %, 13 
as column vectors, 





vol(P) = |det(x1 x2 23), 


the absolute value of the determinant of the 3 x 3 matrix formed from the column vectors. 
This formula has the advantage of being easy to calculate. But it suffers from two related 
disadvantages: (1) it depends on the coordinates of the edge vectors and thus is not very 
geometric, and (2) it requires the ambient space to be R°. 
Both of these objections can be overcome by resorting to the by now familiar device of 
computing the matrix product 





(1 £2 £3)'(1 22 03) = | x8 | (x1 22 23) 


U1, Oet, G1 OG Ge 
20% 20% 
T3ZO0%T1 THOTQ 30 V3 


| 
8 
bo 
e 
8 
e 


Notice that whereas both matrices on the left side of this equation depend on the coordi- 
nates of 71, 2, and x3, the matrix (2; ex;) on the right side depends only on the (geometric) 
inner products of the vectors! The properties of determinant produce the formula 


(det(x x2 #3))? = det(a;e x;). 


Thus we have obtained the formula 


vol(P) = 4/det(x; e x;). 
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This is our desired formula. It displays vol(P) in such a way that we no longer need the 
assumption P C R®. For if the ambient space is RY, we can simply regard 71, v2, 73 as lying 
in a 3-dimensional subspace of R‘ and use the formula we have just derived. The point is 
that the inner products 2,2; from R% are the same inner products as in the subspace of RN, 
since the inner product has a completely geometric description. 





















































e Dimension 2. Here we are going to present the procedure for obtaining the formula 
analogous to the above, but without using a coordinate representation at all. We shall see in 
Section C that this procedure works for all dimensions. 





We suppose P to be the 2-dimensional parallelogram in R% with vertices 0, 71, v2, and 
XL1 +22. We choose the “base” of the parallelogram to be [0,2,| and then use the classical 
formula for the area of P: 


. . x9 
base times altitude 


vol(P) 


Ilara|h h 


x} 


(see the figure). The altitude h is the distance from x2 to the “foot” tx, of the line segment 
from 2x2 orthogonal to the base. The definition of t is then that 


(tz, —2_)e2r, = 0. 
(How many times have we used this technique in this book?!) The definition of h then can be 
rewritten 
h? = |lay—tay|[? = (x2 — tay) @ (29 — tay) 
= (%2—tx1)@ 2p. 


We now regard what we have obtained as two linear equations for the two “unknowns” t 
and h?: 


tz, e2,+0h? =2, 029, 
tt, e%. +h? =220 2». 


We don’t care what t is in this situation, so we use Cramer’s rule to compute only the desired 
me 
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LT, OU, Le XM 
det 
U1, e@% Ge 


det (“1°71 : 
X1e@ XQ 1 
det (a; @ x;) 
Ilr a|/? 


vol(P) = 4/det(x; @ x;). 


This is the desired analog to what we found for n = 3, but we have obtained it without using 
coordinates and without even thinking about R?! 


ie = 





Thus 








EXAMPLE. We find the area of the triangle in R? with vertices 7, 7, k. We handle this 
by finding the area of an associated parallelogram and then dividing by 2 (really, 2!). A 
complication appears because we do not have the origin of R® as a vertex, but we get around 
this by thinking of vectors emanating from one of the vertices (say 7) thought of as the origin. 
Thus we take 





1 0 


= t-7 = —1 ,% =k-7= —l 
0 1 


Then the square of the area of the parallelogram is 


det(a;ex;) = det é 
3. 





Thus, 
1 
area of triangle = 5V3. 














PROBLEM 8-2. Compute the area of the triangle in R3 whose vertices are a2, bj, and 
ck. Express the result in a form that is a symmetric function of a, b, c. 


B. The Gram determinant 
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Our work thus far indicates that the matrix of dot products (2; ex;) will play a major role 
in these sorts of calculations. We therefore pause to consider it as an object of tremendous 
interest in its own right. It is named after the Danish mathematician Jorgen Pedersen Gram. 


DEFINITION. Suppose 21, x2,...,%n are points in the Euclidean space R*. Their Gram 
matriz is the n x n symmetric matrix of their inner products 


G=(xj,erj) = 





UL OT %C1eOeXM~ =... G1 OM 


Ln OL, In@XQ ... IneLn 


The Gram determinant of the same points is 


Gram(21,22,... 


Here are the first three cases: 


iq) = detG 
= det(x; ¢ x;). 


Gram(z) = zex = ||z\|’; 


2 
Gram(x 1,22) = det ( jz" rie a 


z ex, ||xol|? 


= |[e1||*l|z2l|? — («1 ¢ 22)”; 


I|x1 ||? T1OXQ 1X3 


Gram(21,22,%3) =det | z,¢2_ ||x9|\?_ zr. @ rg 


L1 OZ L203 I|r3||? 


Notice that for n = 1, Gram(xz) > 0 and Gram(z) > 0 unless x = 0. For n = 2, 
Gram(21,22) > 0 and Gram(2, 22) > 0 unless x; and 2 are linearly dependent; this is the 
Schwarz inequality. We didn’t encounter Gram(21, 22,23) until the preceding section, but 
we realize that the analogous property is valid. In fact, we shall now give an independent 
verification of this property in general. 


THEOREM. /f xj,.. 


particular, 


0 & 





R, then their Gram matrix G is positive semidefinite. In 


Grama i534 ..qat,) 20. 


6 Chapter 8 


Moreover, G is positive definite <=> 
Gram (#1, :2-;%,) > 0 


<> 11,...,%y are linearly independent. 


PROOF. We first notice the formula for the norm of 4 tiXi: 


| Yan i 


| 
= 
e 
a 
8 
S&S 


i=1 j=l 
n 
ij=1 


In terms of the Gram matrix, the calculation we have just done can thus be written in the 


form : 
| So tia: |P = Gtet, 
i=l 


where ¢ is the column vector with components ¢1,..., tn. 

We see therefore that Gtet > 0 for all tf € R”. That is, the matrix G is positive semidefinite 
(p. 4-26). Thus det G > 0, thanks to the easy analysis on p. 4-26. (Remember: the eigenvalues 
of G are all nonnegative and det G is the product of the eigenvalues.) This proves the first 
statement of the theorem. 

For the second we first note that of course G is positive definite <— > detG > 0. If these 
conditions do not hold, then there exists a nonzero t € IR” such that Gt = 0. Then of course 
also Gt et = 0. According to the formula above, > t;x; has zero norm and is therefore 
the zero vector. We conclude that 21,...,2, are linearly dependent. Conversely, assume 
yo titi = 0 for some t 4 0. Take the dot product of this equation with each x; to conclude 
that )~"_, t;a,e7; =0. This says that the j‘" component of Gt = 0. As 7 is arbitrary, Gt = 0. 
Thus det G = 0. 

















QED 


PROBLEM 8-3. For the record, here is the general principle used above: prove that 
for a general real symmetric n x n matrix A which is positive semidefinite and for t € R”, 
Atet=0 <> At=0. 














The following problem strikingly shows that Gram matrices provide a concrete realization 
of all positive semidefinite matrices. 
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PROBLEM 8-4. Let A be a real symmetric positive semidefinite n x n matrix. Prove 
that A is a Gram matrix; that is, there exist 71,...,2, € R” such that A = (x; e 2;). 
(HINT: use VA.) 





PROBLEM 8-5. Let A be a real symmetric positive definite n x n matrix. Suppose 
A is represented as a Gram matrix in two different ways: 
A = (xa,e2x;) = (ye y;). 


Prove that 
i= By, 10, 


where B is a matrix in O(n). Prove the converse as well. 


As we noted earlier, the inequality Gram(z1, 72) > 0 is exactly the Schwarz inequality. 
We pause to demonstrate geometrical significance of the next case, Gram(21,%2,23) > 0. 
Expanding the determinant, this says that 


ll ||? |]z2|?|l@sll? + 2(x1 © £2) (2 © ©3)(L3 © 21) 


— I|z1||?(z2 . 3)" = I|z2|? (ws ° ii) = I|z3||? (x1 ° 2)” = 0. 





Now consider points which lie on the unit sphere in R‘. These are points with norm 1. For 
any two such points x and y, we can think of the distance between them, measured not from 
the viewpoint of the ambient R™, but rather from the viewpoint of the sphere itself. We shall 
later discuss this situation. The idea is to try to find the curve of shortest length that starts 
at x, ends at y, and that lies entirely in the sphere. We shall actually prove that this curve 
is the great circle passing through x and y. Here’s a side view, where we have intersected R‘ 
with the unique 2-dimensional plane containing 0, x and y (minor technicality: if y = +a then 
many planes will do): 
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The angle @ between x and y is itself the length of the great circle arc joining x and y. This 
is of course calculated by the dot product: cos@ = rey. As —1 < xey <1, we then define 
the intrinsic distance from x to y by the formula 


d(x,y) = arccosxe y. 


(Remember that arccos has values only in the interval [0,7].) The following properties are 
immediate: 


O0<d(z,y) <7; 
2) = SS Sy 
U(x,y) = d(y, x). 
It certainly seems reasonable that the triangle inequality for this distance is true, based 
on the fact that the curve of shortest length joining two points on the sphere should be the 


great circle arc and therefore should have length d(x, y). However, we want to show that the 
desired inequality is intimately related to the theorem. The inequality states that 


d(x,y) < d(x, z) + d(z,y). 
We now prove this inequality. It states that 
arccos Ze y < arccos x e@ z + arccos z @ y. 


If the right side of this equation is greater than 7, then we have nothing to prove. We may 
thus assume it is between 0 and 7. On this range cosine is a strictly decreasing function, so 
we must prove 

xey > cos(arccos x @ z + arccos z @ y). 
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The addition formula for cosine enables a computation of the right side; notice that 


sin(arccosrez) = \/1—(rez)?. 


Thus we must prove 





rey > (vez)(zey) — V1 (weapPyi— (ey. 


That is, 





Vi- (weapy1— (oy > (02)(z0y) —cey. 


If the right side is negative, we have nothing to prove. Thus it suffices to square both sides 
and prove that 


(1 — (we z)’)(1—(zey)”) > ((wez)(zey)—wey)’. 


A little calculation shows this reduces to 


1+ 2x 02)(z0y)(eey) —(we2)?—(z0y)?— (wey)? >0. 


Aha! This is exactly the inequality Gram(z, y, z) > 0. 
PROBLEM 8-6. In the above situation, show that equality holds in 
dist(z,y) = dist(z, z) + dist(z, y) 


if and only if z lies on the “short” great circle arc joining x and y. 
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PROBLEM 8-7. (Euclid) Show that when three rays in R® make a solid angle in 
space, the sum of any two of the angles between them exceeds the third: 


—. 
al? 


PROBLEM 8-8. Suppose m rays meet in convex fashion at a point p in R?. Show 
that the angles at p sum to less than 27. 
(HINT: sum all the angles a, 3, @ over all the triangles. Use the preceding problem.) 
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PROBLEM 8-9. Here is a familiar and wonderful theorem of plane Euclidean 
geometry. Consider the angles shown in a circle x 
with center 0. Both angles subtend the 
same arc of the circle. We say a is 
a central angle, 3 an inscribed angle. y 
The theorem asserts that a = 2/3. 

Now label the three points on the 
circle as shown and show that the equation 
a = 2G is exactly the equation Gram(z, y, z) = 0. 





a] 


C. Volumes in all dimensions 
We are now ready to derive the general formula for volumes of parallelograms. We shall 
proceed inductively, and shall take an intuitive viewpoint: the n-dimensional parallelogram P 
with edges 21,...,2n has as its “base” the (n — 1)-dimensional parallelogram P’ with edges 
X1,.-.,2p—-1. We then compute the altitude h as the distance to the base, and define 
vol,(P) = Voln—1(P’)h. 


Xn 


“Side view”: 


O P’ y 


This “definition” is not really adequate, but the calculations are too interesting to miss. 
Moreover, Fubini’s theorem in the next chapter provides the simple theoretical basis for the 
validity of the current naive approach. 

To complete this analysis we simply follow the outline given in Section A for the case 
n = 2. That is, we want to investigate the “foot” y = eS t;a; of the altitude. We have 

(y—Ip)ex;, = 0 for 1<i<n-1. 
These equations should determine the t;’s. Then the altitude can be computed by the equation 
he = lan — yl? 
(Zn — y) © (Ln — y) 
(Ln — Y) @ Lp. 
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We now display these n equations for our “unknowns” fy,...,t,—1, h?: 
n—-1 
Yo tape xj; +0h? =x, erp, t<2 n= 1: 
j=l 
n—-1 
tjtne rj +h? =Ty_e Ly. 
j=l 


Cramer’s rule should then produce the desired formula for h?: 








Ly, OM, «1. LL OM, VOM 
H2eOXM, ... LOC, G2O My 
det 
12 — LPL ad - yO Teer Bee 
U1, e@U, ... 2102p, O 
TeX, ... LoeLp-1 O 
det 
Ine, ... LneLpn_-1 1 
Aha! This is precisely 
pe Gram(21,...,2n) 
Gram(%1,...,%n—1) 
Of course, this makes sense only if the determinant of the coefficients in the above linear 
system is nonzero. That is, only if 7,...,2%,_1 are linearly independent. Notice then that 
h? = 0 =>} 2,...,2n are linearly dependent; that is, = x, lies in the (n — 1)-dimensional 


subspace containing the “base” of the parallelogram. 
Thanks to this formula, we now see inductively that 














voln(P) = /Gram(z1,...,2n). 





PROBLEM 8-10. In the special case of an n-dimensional parallelogram P with edges 
L1,.--,Xyn € R”, show that 














vol,(P) = |det(#y ay 2.2 @,)|: 


D. Generalization of the cross product 
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This is a good place to give a brief discussion of “cross products” on R”. We argue by 
analogy from the case n = 3. Suppose then that 71,...,2%,_1 are linearly independent vectors 
in R”. Then we want to find a vector z which is orthogonal to all of them. We definitely 
anticipate that z will exist and will be uniquely determined up to a scalar factor. 

The formula we came up with for the case n = 3 can be presented as a formal 3 x 3 
determinant as 














z= det{[ xz, 2 


FDowWo So 


where x; and x2 are column vectors and we expand the “determinant” along the third column. 
This is a slight rearrangement of the formula we derived in Section 7A. 
We now take the straightforward generalization for any dimension n > 2: 


z = det] x, 2... Lp_1 


A 


En 


This is a formal n x n determinant, in which the last column contains the unit coordinate 
vectors in order. The other columns are of course the usual column vectors. 
Notice the special case n = 2: we have just one vector, say 


Xx 
QM 


and the corresponding vector z is given by 


_ a E1 
z = det é >) > x, 


= aey = bey 


are) 








NOTATION. For n > 3 we denote 


Z = 4&1 XXQK +++ KX Mn-1,; 
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and we call z the cross product of the vectors 71,...,%p—,. For n = 2 this terminology doesn’t 
work, so we write 
2 = Ja, 


where J can actually be regarded as the 2 x 2 matrix 


0 -1 
fe a 

(Another way to think of this is to use complex notation x = a+ ib, and then Jx = ix 
(= —b+ia).) If we depict R? as in the above figure, then J is a 90° counterclockwise rotation. 
(Thus J is in SO(2).) 

We now list the properties of the cross product. We shall not explicitly state the corre- 
sponding easy properties for the case n = 2. 
PROPERTIES. 





1. v1 X +++ X Xp_1 is a linear function of each factor. 
2. Interchanging two factors in the cross product results in multiplying the product by —1. 


3. The cross product is uniquely determined by the formula 














(U1 X-++ XAp-1)ey = det(x, %o°--Xn_-1 y), ally e€ R”. 


A. (41 X +++ X Ln-1) @ ZH = O. 





By. ln 3 eo Daca P= det Ge ey at = K ac 

6. 21 X-++ X Lp_-1 =O S > 1},..., Ln_1 are linearly dependent. 
(Proof: < follows from 5, as the right side of 5 is zero. To prove the converse, assume 
X1,...,£n_1 are linearly independent. Choose any y € R” such that 21,...,2%n_1, y are 











linearly independent. Then 3 implies (a1 X...X Zn_i)ey #0. Thus 71 X...X Ln_1 # 0.) 





T.. |e Xo? * Gel] —4/ Grama... 2a-a). 
(Proof: let z= x, X +++ X &p_1. Then 5 implies 


Ilz||* (det(x1...2n_1 2))? 


= (Grail Fi56..5 yas). 
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Now in the computation of this Gram determinant all the terms x; e z = 0. Thus 


U1 OU, wwe LY OL 1 0 
|z\|* = det ; 
In-1@ Ly... Ln-1 @ Xy_4 0 
0 ar 0 Zez 
= |\z||?Gram(a1,...,2n_1).) 
8. If 71,...,@,_1 are linearly independent, then the frame 
{@1,--+,;@n-1,%1 X... X Ln_1} has the standard orientation. 


(Proof: the relevant determinant equals ||z, x ... X @p—1||? > 0, by 5.) 


This completes the geometric description of the cross product. Exactly as on p. 7-9, we 
have the following 





SUMMARY. The cross product 71 X...X %n_1 of n—1 vectors in R” is uniquely determined 
by the geometric description: 


(1) ay % ss. ® My 18 Orthogonal t0.%j 5... 215 

(2) |lay X ... X &p_1|| is the volume of the parallelogram with edges 21,...,%n—1, 

(3) in case 21 X... X Zn_-1 #0, the frame {21,...,%n-1, V1 X ... X Ln_1} has the standard 
orientation. 


PROBLEM 8-11. Conclude immediately that 


é, X... X €, (with é; omitted) = (—1)"*. 


PROBLEM 8-12. Now prove that 


(v1 X...X Ln-1) @ 61 X... X En (with é; omitted) 
= det((%1...%p_1) (with row 7 omitted)). 
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PROBLEM 8-13. Now prove that 


(aX... X Bn_1) © (yi X ... X Yn-a) = det (a; © y;). 


(HINT: each side of this equation is multilinear and alternating in its dependence on the 
y;'s. Therefore it suffices to check the result when the y;’s are themselves unit coordinate 
vectors. ) 





The result of the last problem is a generalization of the seventh property above, and is also 
a generalization of the Lagrange identity given in Problem 7-3. 


PROBLEM 8-14. In the preceding problem we must be assuming n > 3. State and 
prove the correct statement corresponding to R?. 


